QUANTUM PHASE TRANSITIONS IN SPIN SYSTEMS 
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We discuss the influence of strong quantum fluctuations on zero-temperature phase 
transitions in a two-dimensional spin-half Heisenberg system. Using a high-order 
coupled cluster treatment, we study competition of magnetic bonds with and with- 
out frustration. We find that the coupled cluster treatment is able to describe the 
zero-temperature transitions in a qualitatively correct way, even if frustration is 
present and other methods such as quantum Monte Carlo fail. 

1 Introduction 

Phase transitions have been a subject of great interest to physicists over many 
decades. Besides classical or thermal phase transitions, the so-called quantum 
phase transitions (or zero-ptemperature transitions) have started to attract a lot 
of attention (e.g., see Rem). The study of a variety of simple models allows us 
to understand which aspects of thermal and zero-temperature phase transitions 
are common to classes of models and which are more special. For continuous 
order-disorder transitions we basically need the interplay between the interparticle 
interactions and fluctuations. For thermal transitions the Ising model may serve 
as the simplest model. The equilibrium state corresponds to a minimum of the 
free energy, and we have competition between energy and entropy controlled by the 
temperature. For zero-temperature transitions no thermal fluctuations are present, 
and the fluctuations arise due to Heisenberg's uncertainty principle. A correspond- 
ing basic model which has strong quantum fluctuations is the spin-half Heisenberg 
antifcrromagnet (HAFM), particularly in low dimensions. 

The subject of quantum spin-half Heisenberg antiferromagnetism in low- 
dimensional systems has attracted a great deal of interest in connection with the 
magnetic properties of the high-temperature superconductors. Although we know 
from the Mermin- Wagner theoremH that thermal fluctuations are strong enough 
to destroy magnetic long-range order (LRO) at any finite temperature in ID and 
2D, the role of quantum fluctuations is less understood. It is now clear that the 
ground-state of the HAFM in ID is not long-range ordered, whereas the HAFM on 
the square lattice is long-range ordered (e.g., see Ref. 13). However, in 2D there are 
many other lattices with different coordination numbers and topologies, and there 
is no general statement concerning zero-temperature Neel-like LRO. 

Anderson and Fazekas have suggested^ that additional competition between 
magnetic bonds may increase quantum fluctuations and can suppress the Neel-like 
LRO in 2D. Indeed, the strength of this competition may serve as the control pa- 
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rameter of a zero-temperature order-disorder transition. The competition between 
magnetic bonds in quantum spin systems can be caused in various ways. As for 
classical spin systems, frustration can affect the magnetic ordering in quantum spin 
systems. In the classical HAFM the frustration often leads to canted (e.g., spi- 
ral) spin states which may or may not have counterparts in the quantum HAFM. 
Furthermore, due to frustration Marshall's sign rule need not be fulfilled!] The 
violation of the sign rule in frustrated systems makes their theoretical investiga- 
tion particularly difficult. For example, the quantum Monte Carlo (QMC) method 
suffers from the minus sign problem in frustrated spin systems. 

A generic model of a frustrated HAFM is the spin-half J1-J2 model on the square 
lattice, where the frustrating J2 bonds plus quantum fluctuations yield a second- 
order transitiorkfrom a N eel-ordered state to a disordered quantum spin liquid (see, 
e.g., Refs. BlJBb). On the other hand, there are examples where frustration leads 
to a first-order transition in quantum spin systems in contrast- .to .a .second-order 
transition in the corresponding classical model (see, e.g., Refs. t3'tlHl3't3). 

Besides frustration there is a second type of competition between bonds which 
favours a Neel-like distribution of spin correlations over the lattice and other bonds 
which favour the formation of local spin singlets. By contrast to frustration, which 
yields competition in quantum as well as in classical systems, this type of competi- 
tion is present only in quantum systems. The formation of local singlets is accompa- 
nied by the 'melting' of the magnetic LRO. This mechanism for breaking magnetic 
LRO may be relevant for the quantum disordered state in bilayer systemsEJllj as 
well as in CaV^Og (see, e.g., Refs. EjO). Of course, halh-mechanisms can be mixed 
as, for instance, in SrCu 2 (B0 3 )2 (see, e.g., Refs. BBBB). 

In this paper we discuss the quantum order-disorder transition driven by lo- 
cal singlet formation as well as the influence of quantum fluctuations on zero- 
temperature transitions driven by frustration. To that end we study a spin-half 
model on the square lattice in which both mechanisms, frustration and singlet 
formation, are observed in different regions. High-order implementations of the 
coupled cluster method (CCM) are used to obtain a consistent description of both 
types of competition for this model. The CCM (see, e.g., Ref. 123) is one of the 
most powerful and most universal techniques in quantum-^raJaody theory, and 
has previously been applied to quantum spin systemsSEaEafltiaEa with great suc- 
cess. In particular, we shall study to what extent the CCM is able to describe 
zero-temperature transitions in spin systems. 

2 The Model 

We consider a spin-half Heisenberg model on a square lattice with two kinds of 
nearest- neighbour bonds J and J', as shown in Fig. [j], 

H = J^Si-Sj-h J'^Si-Sj. (1) 

The expressions and (2,7)2 indicate nearest-neighbour bonds arranged in a 

regular zigzag pattern, as shown in Fig. [l] by the dotted and solid lines, respectively. 
Each square-lattice plaquette consists of three J bonds and one J' bond. If J' and J 
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Figure 1. Illustration of the model (J-bonds correspond to dotted lines and J'-bonds to solid 
lines) and of the classical spiral state for antiferromagnetic J = +1 and ferromagnetic J' < —1/3 
(left graph) and vice versa (J = — 1, J' > 1/3) (right graph). As discussed in the text, both spiral 
states can be transformed into each other by reversing all of the spins on the B sublattice. 

have different signs then the plaquettes are frustrated, whereas competition without 
frustration is realized for antiferromagnetic bonds J' > J > 0. 

3 The Classical Ground State 

To discuss the influence of quantum fluctuations on the ground-state (GS) prop- 
erties of the model we need to know the classical GS of Eq. (|l|). We set J = 1. 
Without frustration (J' > 0), the classical GS is the Neel state, independent of the 
strength of J'. Frustration appears for (ferromagnetic) J' < 0. One finds that the 
Neel state remains the GS for small amounts of frustration, 0>J'>— 1/3. At 
the critical point J' = —1/3, a second-order transition takes place from the Neel 
state to a spiral state characterized by a pitch angle <& c i = arccos(y / l — l/J'/2) 
(see left graph of Fig. 0). We note that $ c i = (for J' > —1/3) corresponds 
to the Neel state. The spin directions belonging to the A and B sublattices re- 
spectively, are given by Sa (R) = e x cos Q R + e y sin Q • R and S_b (R + x) — 
e x cos(Q • R + 7r 4- 3$) + e y sin(Q ■ R + n + 3$), where R runs over the sites of 
the sublattice A, and Q = (2$,0). We note that this spiral state is incommensu- 
rate in the ^-direction. By contrast to the quantum case, the classical model with 
J = — 1 can be transformed into the model with J = 1 considered above by the 
simultaneous substitution J — > — J, J' — > — J', S ie B — > — Si S s- Hence the physics 
for J = —1 is classically the same as for J = +1 (c.f., Fig.Q). 

4 Competition Without Frustration 

In this section we restrict our attention to the region where J = +1 and J' > 0. 
The terms within the Hamiltonian now compete because the first term favours Neel 
LRO on the "honeycomb" lattice, whereas the second term favours an uncorrelated 
product state of local pair singlets (see below). Again, we emphasize that there is 
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Figure 2. Sublattice magnetization versus J' (for J = 1) (left graph) and energy versus sublattice 
magnetization (right graph) using a mean field approach (g). 



no competition in the classical model. 

Mean Field Approach: We start with a simple mean-held (MF) like description 
of the order-disorder transition. The corresponding uncorrelated MF state for Neel 
LRO is the Neel state |<^MFi) = | TIT ■ ■ •); an d for the dimerized singlet state it is the 
rotationally-invariant product state of local pair singlets \4>mf 2 ) — TI^aH Viii+x) 
~ I liTi+i)]/v2 where i and i+x correspond to those sites which cover the J' bonds. 
In order to describe the transition between both states, we consider anjipcorrelated 
product state interpolating between |</>MFi) and \4>mf 2 ) of the formQo 

I*mp(<)) = II ~n=* [I T<ii+*> - t\ UU+*)] • (2) 
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We have |*mf(* = 0)> = |0 M Fi) and |# M f(* 
(^mfI^I^mf) with respect to t and obtain 



1)) = |0mf 2 )- We minimize 
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(3) 



for the energy per site. For the sublattice magnetization m = (^mfIS^j^^mf) we 
get m = y/ (3 J — J')(3 J + J') / (6 J) for J' < 3 J and m = otherwise. Note that m 
vanishes at a critical point J' c = 3J, and that the critical index is the MF index 1/2. 
Eq. (0) may be rewritten in terms of m as Emf/N = — gJ' — — 4m 2 — 

| J 77? , and Fig. || illustrates that the dependence of Emf on m corresponds to a 
typical scenario of a second-order transition. We can expand Emf up to the fourth 
order in m near the critical point and find a Landau-type expression, given by 
Emf/N = — |j' + | ( J' — 3J) m 2 + m 4 . However, as discussed elsewhere 
for a similar magnetic model for CaA^OgjLZl MF theory probably does not describe 
the critical behaviour correctly. . .. . 

CCM: Let us now apply a high-order CCM approach (for details see Refs. 
to this model. We set the classical collinear Neel state to be the reference state |$). 
We calculate the GS wave function, |\&) = e s |<I>) within the LSUBn approximation 
scheme up to n = 8 and extrapolate to n — > oo. The CCM results for the order 
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Figure 3. Sublattice magnetization (left graph) and excitation gap (right graph) versus J'. 

parameter are shown in the left graph of Fig. || and they are compared to results 
of linear spin wave theory (SWT), exact diagonalization (ED) of N=16,18,20,26,32 
sites, and the MF theory. The CCM is able to describe correctly the order-disorder 
transition, whereas conventional SWT cannot (for more details concerning the SWT 
and ED results see Ref. c3). The critical value predicted by extrapolation of the 
LSUBn results is, however, found to be slightly too large. We may also consider 
the inflection points of m versus J' for the LSUBn approximations. It is assumed 
that the true m( J')-curve will have a negative curvature up to the critical point. 
Thus we might expect that (for increasing n) the inflection point approaches the 
critical point. We find the corresponding inflection points at J' = 3.1 (n=2), 
J' = 3.0 (n=4), J' = 2.9 (n=6) and J' = 2.85 (n=8), indicating a critical value 
J' c somewhere between 2.5 J and 3 J. Notice, that the estimation of 2.5 '^J.LLJ <3 
is consistent with results of series expansions and exact diagonalizations.E3c3 The 
breakdown of Neel LRO due to singlet formation is also accompanied by the opening 
of an excitation gap between the singlet GS and the first triplet excitation. This 
behaviour is well described by the CCM (right graph of Fig. |J) which predicts that 
the gap opens in the range 2 J < J' < 3 J (and notice that the non-zero gap below 
2 J is a result of the limited accuracy of the extrapolation). 

5 Competition with Frustration 

We now consider the frustrated model (where J and J' have different signs). Due 
to the incommensurate classical spiral state the ED technique for finite-size systems 
is less appropriate. The CCM intrinsically considers the limit N — > oo from the 
outset and thus has no problems in dealing with incommensurate states. Hence 
the CCM appears to be particularly suitable to attack the frustrated quantum 
model. We choose the classical state to be our CCM model state, although quantum 
fluctuations may change the pitch angle of the spiral phase. Hence, we determine 
the 'quantum pitch angle' $ by minimizing -Elsubk^) with respect to $. 

Neel versus Spiral: We consider (antiferromagnetic) J = +1 and (ferromagnetic) 
J' < 0. Results for £($) and <&(</') are shown in Fig. |[ The main results are 
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Figure 5. Energy versus quantum pitch angle for LSUB4 (left graph) and quantum pitch angle 
versus .}' (right graph). Note, that <E> = corresponds to the fully polarized ferromagnetic state. 



that: (i) In the quantum case the quantum Neel state remains the GS up to much 
stronger frustration than in the classical case. Indeed, it is generally found for spin 
systems, that quantum fluctuations favour collinear spin structures as opposed to 
noncollinear ones, (ii) The quantum fluctuations change the phase transition from 
second order to first order, (iii) The CCM yields a consistent description of the 
collinear and the spiral phases. 

Ferro versus Spiral: We now consider the model for (ferromagnetic) J = — I 
and (antiferromagnetic) J' > 0. In the classical model we again have a second- 
order transition from a collinear to a spiral state (see section |^) . In the quantum 
model the situation is quite different. Although the collinear antiferromagnetic state 
possesses strong quantum fluctuations, the collinear, fully polarized, ferromagnetic 
state possesses no such quantum fluctuations. The corresponding results for E{&) 
and $(</') are shown in Fig. ||. By contrast to the situation at J = +1, the transition 
from the ferromagnetic to the noncollinear spiral is now of second order, the same 
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Figure 6. On-site magnetic moment versus J' for J = +1 (left graph) and J = —1 (right graph). 



as for the classical model. Furthermore, the classical critical point J'/ J = —1/3 
also holds for the quantum case. The difference between both cases also becomes 
evident when the order parameter (Si) is considered (Fig. ||). For J = +1 there 
is a discontinuity in (Si) at every level of LSUBn approximation. However, the 
extrapolation to n —* oo becomes imprecise close to the phase transition point. We 
cannot therefore decide whether the order parameter vanishes near to the transition 
point. For J = —1 there is a smooth change in (Si) at the critical point. Increasing 
the antifcrromagnetic J' the spiral magnetic order becomes weaker and vanishes at 
/« 1. The underlying reason for that is local singlet formation, as discussed in 
section |], and the continuous vanishing of the spiral order is therefore very similar 
to this second-order transition. However, the strength of J' needed for local singlet 
formation is much smaller due to the assisting effects of frustration (c.f., RefJla). 



6 Conclusions 

In this article we have investigated the zero-temperature phase transitions of a spin- 
half Hciscnberg system on the square lattice. The main results of our treatment 
are: (i) Quantum fluctuations plus competition without frustration are able to de- 
stroy Neel LRO by local singlet formation. This is a pure quantum effect and has 
no classical counterpart. The control parameter is the strength of the competition 
and the breakdown of Neel ordering is accompanied by the opening of a spin gap. 
Standard SWT (even to higher orders) fails to describe this transition, whereas 
the CCM describes both the order parameter and the gap satisfactorily. Since we 
have no frustration, most standard techniques (e.g., QMC) am applicable and a 
quantitative description is possible. As was discussed in Ref.tZI, the critical prop- 
erties seem to correspond to the 3D classical Heisenberg model, (ii) Competition 
due to frustration was found to give more complex magnetic properties. In the 
model considered we have a second-order transition between collinear (antiferro- or 
ferro-magnetic) and noncollincar (spiral) states driven by frustration in the clas- 
sical case. In the quantum spin-half model, standard techniques (e.g., QMC) are 
not applicable due to the violation of Marshall's sign rule. By contrast, the CCM 
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provides a consistent description of collinear, noncollinear, and disordered phases. 
Furthermore, we find a strong influence of quantum fluctuations on the nature of the 
collinear-noncollinear transition, and quantum fluctuations (which favour collinear 
ordering) may change the second-order classical transition to a first-order quan- 
tum transition. If quantum fluctuations are suppressed in the collinear phase, the 
transition to the spiral phase is similar for the quantum and classical models. 
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